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Abstract. I describe analysis of correlations in the outputs of the three LIGO 
interferometers from LIGO's first science run, held over 17 days in August and 
September of 2002, and the resulting upper limit set on a stochastic background 
of gravitational waves. By searching for cross-correlations between the LIGO 
detectors in Livingston, LA and Hanford, WA, we are able to set a 90% confidence 



PACS numbers: 04.80.Nn, 07.05.Kf, 98.70.Vc 
E-mail: jtwhelanaioyno.edu 



t Published as Class. Qu antum Gray. 2 1 (2004) S685-S690; 

available online at stacks.iop.org/CQG/21/S685 copyright (C) 2004 IOP Publishing Ltd. 

§ The members of the LIGO Scientific Collaboration are listed in the paper by Allen and Woan for 

the LIGO Scientific Collaboration, in these proceedings: stacks.iop.org/CQG/21/S671 



First LIGO Science Results: Stochastic Background Upper Limits 



2 



1. Introduction 



The LIGO interferometric gravitational wave (GW) detector held its first science run 
(SI) in 2002, from 23 August to 9 September 0. LIGO consists of an interferometer 
(IFO) with 4 km arms in Livingston, LA, USA (the LIGO Livingston Observatory, or 
LLO), called LI for short, and two IFOs, with arms of 4 km and 2 km, in Hanford, 
WA, USA (the LIGO Hanford Observatory, or LHO), called HI and H2, respectively. 
The data were analysed to search for GW bursts |2] , signals from inspiralling neutron 
star binaries periodic signals from a rotating neutron star 0], and stochastic 
backgrounds |5J. This paper summarizes the analysis method and results of the search 
for a stochastic background of gravitational waves (SBGW), which are explained in 
more detail in 



2. Fundamentals of Analysis Method 

A SBGW is assumed for simplicity to be isotropic, unpolarized, Gaussian, and 
stationary. Subject to these assumptions, the SBGW is completely described by its 
power spectrum. It is conventional to express this spectrum in terms of the GW 
contribution to the cosmological parameter f2 = p/p C rit : 

, ,s 1 dp GW f dp GW 

"GW J = -Jj J = TT~ ■ V 1 ) 

p cr it din/ p crit dj 

Note that ^Gw(f) has been constructed to be dimensionlcss, and represents the 
contribution to the overall Oqw per logarithmic frequency interval. In particular, 
it is not equivalent to d£lcw/df ■ Note also that since the critical density p C rit, which 
is used in the normalization of Ogw(/)j is proportional to the square of the Hubble 
constant Hq [B], it is convenient to work with /!j 00 !Jgw(/)i which is independent of 
the observationally determined value of /iioo = ioo km/s/ m p c ' ^ 

The standard method to search for such a background is to cross-correlate the 
outputs of two GW detectors [Tj- If the noise in the two detectors is uncorrelated, the 
only non-zero contribution to the average cross-correlation (CC) will come from the 
SBGW. In the optimally-filtered CC method (described in more detail in [51 151 110|L 
one calculates a CC statistic 

Y= [ dt 1 dt 2 h 1 (t 1 )Q(t 1 -t 2 )h 2 (t 2 ) = [ dfh${f)Q(f)h 2 (f) (2) 



where h\ t2 {t) are the data streams from the two detectors, hi^if) are their Fourier 
transforms, and Q{t\ —t 2 ) (with Fourier transform Q(f)) is a suitably-chosen optimal 
filter. The choice which optimizes signal-to-noise ratio for a constant- f2ow(/) 
background is ||] 

The normalization of the optimal filter is conventionally chosen so that in the presence 
of a SBGW of strength Qgw(/) = ^o, the expected mean value of the CC statistic is 

(Y) = hj oa n T (4) 

| Although hioo is now much more accurately known than it once was, we still work in terms of 
^lOO^GwC/) t° facilitate comparisons with prior results. 



First LIGO Science Results: Stochastic Background Upper Limits 



3 



where T is the duration of the analysed data sets. The expected variance of the CC 
statistic is 

< cm = jj df Pl (f) \Q(f)\ 2 P 2 (f) cc (/ * Pi [ fff (/) ) 1 ■ (5) 
The method is sensitive to backgrounds on the order of 

The sensitivity of this method improves with time and is limited by the power spectral 
densities Pi^(f) of the noise in the two detectors. The factor 

l{f) = d lab df A e ^/n Ax/c pab {n) (?) 

in the numerator of the integral is the overlap reduction function which describes 
the observing geometry. Here P^d{ n ) is a projector onto symmetric traceless tensors 
transverse to a direction n and d° b 2 are the detector response tensors for the two 
detectors. These are the tensors with which the metric perturbation h ao at the detector 
should be contracted to obtain the GW strain h = d ab h a b- If u a and v a are unit vectors 
pointing in the directions of an IFO's two arms, its response tensor is 

dab = -^{u a U b - V a V b ) . (8) 

The overlap reduction function is equal to unity for the case of a pair of IFOs at the 
same location with their arms aligned, and is suppressed as the detectors are rotated 
out of alignment or separated from one another. It also oscillates with frequency 
as correlations are suppressed for detectors whose separation is comparable to or 
greater than the corresponding GW wavelength. Figure^shows the overlap reduction 
functions for combinations of detectors which were operational during SI. 



3. Prior Results 



The previous best upper limit on a SBGW from direct observation with GW detectors 
was /if 00 f2Gw(900 Hz) < 60 [12], set by correlating the resonant bar detectors Explorer 
(in Geneva, Switzerland) and Nautilus (near Rome, Italy). A broad-band limit of 
^ioo^gw(/) < 3 x 10 5 was set using a pair of 'prototype' IFOs |T5|. 

More stringent upper limits can be set on astrophysical grounds. They are 
detailed elsewhere jS] [S] E] > but we mention the bound from big-bang nucleosynthesis 
HOE], which states that a cosmological SBGW is limited by 

r % hlooSlowif) < 10- 5 • (9) 

This broad-band limit implies that any cosmologically interesting SBGW very likely 
lies several orders of magnitude below the existing limits. 



4. Details of Analysis Method 

Each of the three combinations of detectors (Hl-Ll, H2-L1, and H1-H2) was analysed 
separately for CCs. Since the power spectra P\^{f) varied over the course of the 
SI run, the coincident data for each pair of IFOs were divided into 15 min blocks, 
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Figure 1. The overlap reduction function -f(f) for combinations of the LIGO 
Livingston Observatory (LLO) and LIGO Hanford Observatory (LHO) with each 
other and with the GEO600 site. (GEO600 was also operational during SI, but has 
not been included in this analysis because it was considerably less sensitive than 
the LIGO detectors.) Note that the overlap reduction function for correlations 
between the two detectors at LHO (HI and H2) is identically equal to unity. The 
solid lines (at / = 40 Hz and 314 Hz) show the range of frequencies used in our 
analysis. 



and an optimal filter constructed for each such block. To maximize overall signal- 
to-noise ratio [H], we combined the CC statistics from the different blocks using a 
weighting factor of cr t ~^ or , where ctheor is the theoretical standard deviation defined 
in Note that this can be calculated from the individual power spectra, without 
cross-correlating the data. To avoid problems from noisy and presumably less Gaussian 
data, we discarded the 15 min blocks with the highest er t h oor values, corresponding to 
a 5% total contribution to ^ "theor- 

For each block, an optimal filter was constructed with a frequency resolution 
of 0.25 Hz according to the discrete frequency-domain analogue of 10. The range 
of frequencies included in the calculation of the CC statistic was chosen to be 40- 
314 Hz for Hl-Ll and H2-L1, and 40-300 Hz for H1-H2. Given the power spectra of 
the instruments and the expected spectrum of correlations associated with a constant 
!Jgw(/), frequencies outside that range were not expected to improve the sensitivity 
appreciably. Additionally, individual frequency bins associated with cross-correlated 
instrumental noise were omitted from the sum over frequencies (which means the 

^ The power spectra were constructed using Welch's method, with 449 overlapping Hann-windowed 
periodograms averaged to produce a power spectrum estimate with a resolution of 0.25 Hz. 
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optimal filter was effectively set to zero there). These were integer multiples of 16 Hz 
and 60 Hz, as well as a few frequencies which had a coherence over the entire run 
above 0.2, namely 250 Hz for L1-H2, and 168.25 Hz and 168.5 Hz for H1-H2. 

Each 15 min block was divided into ten 90 s segments using a Tukey window which 
consisted of half s Hann transitions on either side of an 89 s flat top. The CC statistic 
was calculated for each segment, using the discrete analogue of the frequency-domain 
form of (J2J, and these were combined to give a CC statistic for the entire block. In 
this way, we were able to obtain, through the measured standard deviation of the ten 
CC statistics within a block, a statistical measure of the error associated with the 
CC statistic for the block. We also estimated the systematic error associated with the 
change in sensitivity and calibration over the course of each block. + We then combined 
all three, appropriately weighted over the whole run, to obtain a total estimated error 
CTtot associated with the point estimate h1 00 Qo calculated from the weighted average 
of all the CC statistics using QJ. In the absence of cross-correlated noise, the 90% 
confidence level upper limit on f^o, the constant value of ^gw(/)) is 

h 2 10Q n Q < h 2 W0 h Q + 1.28?tot • (10) 

5. Results 

The results for the three IFO pairs are summarized in Table ^ There is a 



IFO Pair 


obs time (h:min) 


h 100 Q 


<?tot 


90% CL UL 


H2-L1 


51:15 


0.2 


18 


23 


Hl-Ll 


64:00 


32 


18 


55 


H1-H2 


100:15 


-8.3 


0.9 


N/A 



Table 1. Summary of the point estimate h^ 00 Qo and total estimated error ertot 
for the three IFO pairs considered. Note that no upper limit is set from H1-H2, 
the two IFOs at the Hanford site, since there was evidence of cross-correlated 
noise. Associated with each of the values quoted is an additional 20% uncertainty 
arising from the calibration of the instruments. 

statistically significant anti-correlation observed between HI and H2, two IFOs which 
share the same vacuum envelope at the LHO site in Hanford, WA. Time-shift and 
X 2 analyses show that this anti-correlation is inconsistent with constant-Slow/ (/) 
SBGW, so we conclude it is due to instrumental cross-correlations between the two 
colocated detectors.* For the inter-site measurements (Hl-Ll and H2-L1), the lack of 
statistically significant cross-correlations makes these checks trivial, and we proceed 
to setting an upper limit from each pair.ft The stronger upper limit is set by H2 and 
LI, and it is 

hlaoClo < 23 ±4.6 . (11) 

+ The calibration of the LIGO detectors was monitored by tracking the amplitude of a sinusoidal 
'calibration line' signal injected into the IFO arm length. The output amplitude, recorded once 
per minute, allows construction of a frequency-dependent response function which accounts for IFO 
alignment drifts, as detailed in [T1I15I . 

* Of course, the fact that it's an anti-correlation rather than a correlation is another reason it can't 
be due to a SBGW. 

(I One might try to combine the different measurements into a single limit, but we choose to consider 
them individually, especially given the observed correlations between the two Hanford IFOs, which 
complicate the issue of combining the Hl-Ll and H2-L1 as supposedly independent measurements. 
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This represents a factor of 2-3 improvement over the previous direct upper limits 
described in Sec. |3| and an improvement by a factor of over 1000 over the previous 
measurements with interferometric detectors. 

For more details on the analysis, the reader is directed to [S]. 
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